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Abstract 

We consider radiative backreaction for global strings using the Kalb-Ramond formalism. 
In analogy to the point electron in classical electrodynamics, we show how local radiative 
corrections to the equations of motion allow one to remove the divergence in the self field 
and calculate a first order approximation to the radiation backreaction force. The effects 
of this backreaction force are studied numerically by resubstituting the equations of mo- 
tion to suppress exponentially growing solutions. By direct comparison with numerical 
field theory simulations and analytic radiation calculations we establish that the 'local 
backreaction approximation' provides a satisfactory quantitative description of radiative 
damping for a wide variety of string configurations. Finally, we discuss the relevance of 
this work to the evolution of a network of global strings and their possible cosmological 
consequences. These methods can also be applied to describe the effects of gravita- 
tional radiation backreaction on local strings, electromagnetic radiation backreaction on 
superconducting strings and other forms of string radiative backreaction. 



1. Introduction 

A variety of unified field theories predict the formation of a network of topological de- 
fects at one or more phase transitions in the early universe [1]. Strings associated with 
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the breaking of local symmetries have generated the most interest in the literature be- 
cause, amongst other reasons, GUT-scale strings could have been the initial seeds for the 
formation of large-scale structure [2]. However, local strings are tightly constrained by 
their contribution to the gravitational radiation background [3]. There are other types 
of strings which circumvent this constraint and which may have similar cosmological 
implications, in particular those formed when a global symmetry is broken. Instead of 
radiating gravitationally, the dominant radiation mechanism for these strings is the emis- 
sion of massless Nambu-Goldstone bosons [4]. In a recent publication [5], we studied the 
nature of this radiation in detail, using analytic and numerical techniques. We demon- 
strated that a low energy effective action known as the Kalb-Ramond action, provided 
an accurate description of the dynamics of global strings even at the moderately high 
velocities one expects in a realistic string network. Within this formalism the topologi- 
cal coupling of the massless field to the string is linearized. One finds that the coupling 
between the field and the string worldsheet is similar to that of the point electron in 
electromagnetism. However, there are still difficulties associated with this approach, 
notably because equations of motion are inconsistent due to a divergent self-field. 

This type of problem has been well understood for some time in the context of a 
point electron in classical electrodynamics [6]. In the case of the electron, the self and 
radiation fields can be distinguished easily since, at large distances R, the self field falls 
off as 1/-R 2 , whereas the radiation field falls off as 1/R. Careful analysis of the equations 
of motion leads to the renormalisation of the electron mass by the Coulomb self-field, 
using the classical electron radius to cut off short distance divergences and a first order 
approximation to the radiation backreaction force known as the Abraham-Lorentz force, 



where X^(r) is the position on the electron's worldline at time r. The dependence of this 
force on has lead to problems in numerical applications since there exist exponentially 
increasing solutions to the equations of motion. These unphysical 'runaway' solutions 
can only be suppressed by rewriting the equations of motion as an integro-differential 




(1) 



2 



equation. 

In a recent letter [7], we proposed a formalism for removing the self force and cal- 
culating a first approximation for the radiation backreaction force of strings using the 
analogy of classical electrodynamics. However, the analogy is not exact because the 
strings are line-like objects of possibly infinite extent. We circumvent this problem by 
assuming that the dominant contribution to the backreaction force comes from string 
segments in the vicinity of the point in question, henceforth known as the 'local backre- 
action approximation'. This approximation will not be valid in every situation, but in 
a wide variety of circumstances it should work well. In this paper we will consider this 
approximation in greater detail. We will elaborate on the derivation of the radiation 
backreaction force and give a discussion of the physical aspects of the approximation. 
Further numerical evidence will be presented in support of the validity of this approx- 
imation in physically important cases by direct comparison between modified Nambu 
dynamics, evolved numerically using the backreaction force, and numerical field theory 
simulations. The one free parameter in our analysis, effectively the damping coefficient, 
can be normalized by comparing with numerical field theory simulations and known 
analytic results. 

One of the main motivations for this work is to implement appropriate radiative cor- 
rections in a full network simulation. We anticipate that the scaling assumption for gauge 
strings, numerically verified in refs. [8,9], will also be seen to be valid for global strings. 
However, it is anticipated that the parameters quantifying the small scale features will be 
somewhat different [10]. Accurate numerical simulations will allow estimates of the cos- 
mological axion density to be refined. A similar formalism is applicable to gravitational 
radiation backreaction on local strings [11] and electromagnetic radiation backreaction 
on superconducting strings. 

Throughout this paper we employ a (H ) signature for the spacetime metric 

g^ u and (H — ) for the induced metric on the string worldsheet 7 a &, the coordinates for 
which are given by = X M (cr, r), with the null coordinates, u = a — r, v = a + r. 
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2. Analytic formalism 



2.1 The Kalb-Ramond action 

Th essential features of global strings in flat space are exhibited in the simple U(l) 
Goldstone model, with action given by 

1 



/ 



(fx 



2\2 



d^d^ - - ft) 



(2) 



where $ = <f>e is a complex scalar field which can be split into a massive (real) com- 
ponent <p and a massless (real periodic) Goldstone boson ft. The analytic treatment of 
global string dynamics is hampered by the topological coupling of the self field of the 
string to the Goldstone boson radiation field. However, we can exploit the well-known 
duality between a massless scalar field and a two-index antisymmetric tensor B^ v to 
replace the Goldstone boson ft in (2) via the relation 

<j> 2 d»ft=\Ue^ Xp d v B x f>. (3) 

Performing this transformation carefully and integrating over the massive degrees of 
freedom about the two-dimensional string worldsheet X M (cr, r) [12], yields the flat-space 
Kalb-Ramond action [13,14], 



S = -hq J dadr^ + ^ J d 4 xH 2 - 2yr/ a J B^da^ , 



(4) 



where H^p = d^B^p + dpB^ +<9 a -Bg M is the field strength of the antisymmetric tensor 
field B^ v , the metric induced on the world sheet is 

1 ab = 9„»d a X»d b X\ 7 = det( 7a6 ), (5) 

and the area element on the worldsheet is 

da^ = e ah d a X»d h X v dadT . 



The first term is the familar Nambu action for local strings, the second is the anti- 
symmetric field strength for both external fields and the self-field of the string and the 
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last term is a contact interaction between the antisymmetric tensor field and the string 
worldsheet. The coupling between the string and the antisymmetric tensor is analogous 
to the electromagnetic coupling of the point electron to the electromagnetic field. This 
analogy underpins our subsequent development of global string dynamics based on (4). 

Varying the action (4) with respect to the worldsheet coordinates and the antisym- 
metric tensor yields the string equations of motion and the tensor field equations, 

r (6) 

diiH i«*f> = -4nJ af3 = -2nf a / dadr5 4 (x - X(a, r))V ap , 

where V aj 3 = e ab d a X a dbXp is the antisymmetric vertex operator. In the conformal string 
gauge and the Lorentz antisymmetric tensor gauge, 

X 2 + X' 2 = 0, X-X' = 0, d^ u = 0, (7) 

the equations of motion (6) become 

UB a p = -4yrJ a/3 = -2yr/ a J dadr5 4 (x - X(a, r))V a p , 



(8) 



where □ = g^d^d^ and V a p = XaX'^—X^Xp. These equations are problematic because 
the self-field diverges as any point of the string is approached, that is x — > X(a, r). 

2.2 Simple string configurations 

If one ignores the effects of the force density then the equations of motion reduce to 
the well-known Nambu equations of motion, a massless wave equation. The equations 
have solution 

X° = t = r, X = i [a(«) + b(v)] , (9) 

where the functions a(it) and b(v) are the left- and right-moving parts of the solution 
(recall u = a — t and v = a + t). Using the conformal gauge conditions (7), one can 
deduce that 

a' 2 = l, b /2 = l. (10) 
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Figure 1: Schematic of the solutions of the Nambu equations of motion: (a) loop solutions 
parametrized by the invariant length L and (b) long string solutions parameterized by the wave- 
length L and the relative amplitude S. 

The equations (9) and (10) have closed loop and periodic long (or infinite) string 
solutions. The loop solutions are parametrized by the length of the loop L, which 
is closely related to the characteristic frequency O = 2tv/L, whereas the long periodic 
solutions are parametrized by the wavelength L and the ratio of amplitude to wavelength 
or the relative amplitude S = 2nA/L, where A is the amplitude. Fig. 1 shows a schematic 
of the two types of solution we shall consider. In general situations, such solutions will 
correspond to a superposition of a large number of harmonics. 

A simple two parameter family of loops, known as Kibble- Turok loops, involve just 
the first and third harmonics [15]: 



where O = 2tt/L and < a < 1, — 7r < ip < n. If a = and ip = then the solution 
is a circular loop, which oscillates between a circle of radius L/2n and a point. For a 
significant range of the parameters a and ip these solutions can be shown to be non 




(1 — a) cos fi-u — -|a cos 30w — cos ip cos Ov, 



(11) 
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self-intersecting and so the dominant decay mechanism is likely to be through radiation 
rather than loop fragmentation. The time evolution of a particular solution with = it/ 3 
and a = 0.5 is shown in fig. 2a. Kibble-Turok loops also generate cusps, that is, points 
on the string which reach the velocity of light. For example, cusps will appear on the 
a = solution at a = L/4, 3L/4 when t = (n + l/2)L/2 (n integer). Whether or not 
cusps are generic on realistic loops has been the subject of various heuristic discussion, 
which have also considered the unknown effect of backreaction on cusp evolution. 

Strings reconnect or 'exchange partners' when they intersect. This process intro- 
duces kinks — contact discontinuities in the velocity X and tangent vector X' — which 
propagate along the string at the speed of light. Because realistic loops are produced by 
long string reconnections or self-intersections we can be sure they will possess at least 
two kinks, probably more. An idealized loop with four kinks between four straight string 
segments is given by the following [16] 

X = 1 [a(* - r) + b(<r + r)] , 

where 

(^>- r )-7) A 0<a-r<n 
a(ff-r) = < ) \ (12) 

(t-^- t) ) a n<a ~ T<2n > 

M* + r)=<)£ iJ (13) 

with arbitrary unit vectors A and B. The two pairs of kinks propagate in opposite 
directions around the loop. In the special case A • B, the loop is planar and oscillates 
between a square and a doubled line. 

A simple, symmetric long string solution can be constructed from equal and oppo- 
sitely propagating helicoidal waves in the fundamental mode [17], 

X = I zrrz [ cos Qu + cos Qv] , — [ sin Qu + sin Qv] , - a/1 — S 2 (u + v) J 

J (14) 



sin Oct cos fit, ^ cos fia cos fit, \J\ — £ 2 o^j , 



where < £ < 1 and £ — > 1 in the relativistic limit. This corresponds to a helicoidal 
solution which oscillates between a static helix and a straight line, as shown in fig. 2(b) 
for £ = 0.6. Because of its perfect symmetry, calculations for the radiation power from 
this solution are analytically tractable [17,5], though such configurations are unlikely to 
be found in a realistic string network. 

One can generalise the helicoidal solution (14) to have unequal left and right moving 
amplitudes, 

X= {^ COSnU+ ^} COsClU ' 2^ Sin ° M+ 2% Sin ° M ' ~ £ R u+ ^\/ 1 ~ £ L v ^j > ( 15 ) 

where < £r,£l < 1- This type of solution is thought to be a reasonably accurate 
description of long strings in a realistic network, since within a sufficiently small volume 
the number of left and right movers are unlikely to be strongly correlated. One special 
case of this solution is that with no left moving perturbation, 



X 



It has been suggested that pure left- or right-moving configurations do not radiate and 
so they will propagate indefinitely [18]. However, we shall argue that such solutions are 
not physically relevant because they require initial data with the string fields artificially 
correlated out to infinity. 

A solution, similar to the helix (14), but with sinusoidal perturbations in only one 
plane is 

X = ( — [ cos fiw + cos Civ] , 0, ^- [E(£, flu) + E(£, flv)] J , (17) 
where E(k, (j>) is the incomplete elliptic integral of the second kind, defined by 

E(k,(f>)= / deVl - k 2 sm 2 6, (18) 
Jo 

where < £ < 1 and £ — > 1 in the relativistic limit. 

Long string solutions of (10) can also possess kinks as in (12). A simple solution for 
a periodic distribution of kinks on a perturbed straight string, consists of the left- and 



right-moving perturbations, X\j and Xy respectively, 

1 



X 



A£ 2 

Xu + X v ,0,^\/1- —^-(u + v) 



(19) 



where 
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(20) 



(21) 



where < £ < ^ and, in this case, £ — > ^ is the relativistic limit. In fig. 2d, notice 
how the kinks split into two which propagate at the speed of light in opposite directions 
along the string. 

2.3 Green functions 

One of the most basic techniques of mathematical physics is the inversion of differential 
equations such as the field equation in (8) using Green functions. The basic Green 
functions satisfy 

UD{x) = 5\x) , (22) 

which implies that the solution to 



□F(x) = S(x) , 



(23) 



is given by integrating the product of the Green function and the forcing term S(x) over 
spacetime 

F(x) = [ d 4 x'D(x - x')S(x') . (24) 



In order to deduce a specific form for the Green function one must specify some 
boundary conditions which define the region over which initial data is known. The two 
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Figure 2: The evolution of various string trajectories: (a) Kibble- Turok loop (11) (ip = 7r/3, 
Ol = 0.5), (b) four kink loop solution (12) (c) helicoidal long string (14) with (E = 0.6), and 
(d) kink solution (19) with 45° openning angle. 

most common Green functions used are the retarded and advanced time Green functions, 
which use initial data on the backward or forward light cones respectively, 

D ret (x) = ±-9(x°)5(x 2 ) , D adv (x) = ^9(-x°)5(x 2 ) , (25) 

where x^ = (x°, x) and 6(x°) is the Heaviside function, that is 6(x°) = 1 for x° > and 
0(x°) = otherwise. 

In problems where radiation is involved one wishes to separate radiative effects 
from those of the self-field. The radiation field is free and the radiation Green function 
must satisfy a homogeneous version of (22). One can construct such a Green function by 
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subtracting the advanced Green function from the retarded. Similarly, one can construct 
the Green function for the self-field by summing the retarded and advanced Green's 
functions. Using appropriate normalisations, one can deduce that 

£>rad = \ (Dret ~ A,dv) = ^e(x°)S(x 2 ) , 

1/ \ 1 (26) 

£>self = g ( A-et + D adv j = ^ S ( x ) , 

where e(x°) = 6(x°) — 6(—x°). One can calculate the self- and radiation-field for a 
problem such as (23) using similar expression to (24) with the basic Green function 
replaced by the appropriate expression from (26). 

2.4 Lienard-Wiechart Potentials 

Using the Green function techniques described in the previous section, one can deduce 
that 

B a p(x) = -Air / d 4 x'D ret (x - x')J a p(x') 

3 , (27) 
= -2nf a / dadr D iet (x - X(a, r))V a p(a, r) . 

The integration is over all time and over all string segments. In the case of a closed loop 

this is a finite range, but for long (infinite) strings the range is infinite. If one defines 

A M = x^ — X^a, f) while treating a as a = <r(f), then 

d(A 2 ) = -2A -Xdf , 



p p d(A 2 ) A- Xdf' 
Substituting into (27) and evaluating the delta function one can deduce the Lienard 

Wiechart potential [6,19,20] 



T=T R 

.21 



(28) 

|A.X| ' 

where A 2 \f= TR = and tr < t. The modulus sign in (28) preserves the orientation of 
the region of integration when evaluating the delta-function. In order to calculate the 
radiation backreaction force one requires the derivative of (28). This can be calculated 
by performing an integration by parts, 



" 1 2 i A.Xdf V A.X J 



(29) 

r=r R 
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One can separate the radiation field from the self field by using the Green functions 
for the self- and radiation- fields (26). Therefore, one can calculate the Lienard-Wiechart 
potentials and their derivatives for both the self and radiation fields. 



Blf(x) = J da 



d,B° a f(x) 



fa, 

T 




A-Xdr \\A-X\ 

1 d ( AfjVcfP 

A-Xdf \\A-X\ 



+ 



1 a (A„y r 



T = TR 



T = T R 



A-Xdr \\A-X\ 

1 d ( AjjYg£ 

A-XdT VIA -XI 



(30) 



where A 2 \ f= 



and tr < t, r' R > t. Effectively, then, we have performed the 



separation of the self- and radiation-fields. An attempt was made to perform this split 
in ref. [21], using techniques similar to those used in classical electrodynamics [6]. This 
method performed the split on the basis of asymptotic fall off. As already discussed, this 
procedure works in the case of the electron since the self- field falls off like 1/-R 2 , whereas 
the radiation field falls off like 1/R for large R. However, this procedure may work for 
string loops. However, it is doomed to failure for long strings, since both the self- and 
radiation-fields fall off as 1/R. 

2.5 The 'local backreaction approximation' 

It has already been noted that the renormalisation procedure for strings is more com- 
plicated than that for the point electron. The main problem becomes obvious when one 
compares the Lienard-Wiechart potentials for strings to those for the electron [6] . Since 
the string is an extended object, the Lienard-Wiechart potential is an integral along the 
string. In the case of a loop of length L, this integral will be in the range < \a — a\ < L 
for a point X{a, r) on the string and can be easily approximated. However, in the case 
of a long (infinite) string the range is — oo < \a — a\ < oo and the integral cannot be 
evaluated without the solution being periodic. 

In more general situations this is not possible and one must make what we shall 
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call the 'local backreaction approximation'. Since the effects of backreaction from string 
segments at large distances from the point in question must be suppressed, it seems 
sensible to truncate the integrals of (28) and (30) at some renormalisation scale A, which 
is at present arbitrary. That is the integrals are over the range —A/2 < \a — a\ < A/2. 
Using the case of the loop of length L as an example seems to suggest that A ~ L. In 
fact our expectation is that, for more general string trajectories, an appropriate choice 
for A would be near the average curvature radius of the string. 

Using the local backreaction approximation one can perform the renormalisation of 
the self-field and the derivation of a first order approximation to the radiation-field. If 
one allows x — > X(a, r), then (30) can be expanded in terms of s = a — a and t = r — f . 
This procedure requires that the natural scale for the otherwise arbitrary renormalisation 
cut-off A be less than the average curvature radius of the string. In this case one finds 
that the condition A 2 | TH = 0, implies that t = \s\ + 0(s 4 ) and A 2 | T ^ = 0, implies that 
t = —\s\ + 0(s 4 ). Ignoring terms of order four in s and t, allows one to deduce that 

/a - \x [f ,V af3] - X'^V aP \ log(A/5) + 0(A 2 ) , 

(31) 



rrself 



trrad 



2X 4 
fa 

~2X 4 



1 • 



(2 X X \ 



X hl V r 



A + 0(A 2 ), 



where ^4[ Ma /3] = + Ap^ a + A a p^. Note that the self-field has no order A term. 

Ignoring terms of order A 2 , we can then obtain expressions for the self- force and the 
first order approximation to the radiation backreaction force density, 



Ff' = -2nfl\og{A/5) 



+ 



AfX.X\ JX.X 
3V^^/ + \~X^ 



+ 2 



X'.X 
X 2 



x u + 



4 fX 
3l J2 



(32) 



where 5 (<< A) is the width of the string core and corresponds to the ultra-violet 
renormalisation scale. These expressions for the self and radiation force densities are 
extremely complicated, however, our confidence that these are the correct expression is 
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strengthened since they non-trivially respect the conformal gauge conditions, that is, 
jrseif . x = 0, J^ self • X' = 0, F™ d • X = and F*** ■ X' = 0. 

The self-field is a multiple of the left hand side of the equations of motion and facil- 
itates the well-known renormalisation of the string tension, in a way exactly analogous 
to the mass of a point electron. The equations of motion in this case are 



MA) 



, (33) 



where /z(A) = + 27r/f log(A/<5) is the renormalised string tension 

For general string trajectories, similar to those discussed in ref. [5] , some of the terms 
in (32) can be shown to be sub-dominant. In particular, if the string solution is specified 
by the relative amplitude £ and its wavelength L, then one finds that X ~ 0{EL) and 
each subsequent derivative requires a division by L. Most of the higher order terms in £ 
can be dropped, though it is necessary to keep two of the higher order terms to maintain 
the gauge conditions. One then finds that it is possible to approximate and F™ A 

by 



A* 



3X 4 

rad „ 47T/ a 2 A 



2 M ... /, «V , , vl , ■ 04) 

- X» - I 

1 -X 2 , 



2.6 Generalization to the temporal transverse gauge 

For flat-space string dynamics, the conformal string gauge is usually employed. However, 
when considering problems in which the string energy decays, it more convenient to use 
the temporal transverse gauge in which X° = t = r and X.X' = with X^ = (£, X). In 
this gauge, the equations of motion for the string (8) are 

^X-i^))=f, fi e = f°, (35) 

where e 2 = X 2 /(l — X' 2 ) (not to be confused with the relative amplitude S) and = 
(/°,ef + /°X). Radiative damping will naturally be incorporated in the decay of the 
coordinate energy density e, rather than in the non-intiutive time redefinitions of the 
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conformal gauge. The string energy and momentum per unit length in the temporal 
transverse gauge are then given by 



E=^ dae, 
L Jo 



f 1 

' L J 



cfcreX . 



(36) 



The renormalisation procedure in this gauge is similar to that for the conformal 
gauge, but with the added complication that X 2 + X' 2 ^ 0. After a detailed set of 
manipulations one can deduce that 



rrself _ /i 



H 



rad 
Ata/3 



2X 4 

fa 



1 - 1 „ 



log(A/<5) + 0(A 2 ) , 



2X 4 



- ^X^V aP] - -X [fl V aP] + —X'^Vap] 



(37) 



'2X.X T/ 



A + 0(A 2 ), 



where ^4[ Ma /3] = + A piia + A aPll . Note that, once again, the self- field has no order 

A term. Ignoring terms of order A 2 , we can then deduce expressions for the self- force 
and the radiation backreaction force, 



f selt = - 2nf 2 log( A/5) ( X - i 



pad =7r/ , A< j _ eX + 



2e 



XX 

1 -X 2 



+ 3e 



e \ 1 - X 5 



+ 



4 / X' • X 



3eVl-X 2 7 e (1-X 2 ) 2 
/0) self = _ 27r / a 2 log(A/5)6, 



4(X-X)(X'-X) e/'X'-X^ 8/ / X X' 



e 2 VI 



X 2 



e 4 VI 



X 2 



/0) rad =7r/ 2 A< -r 



XX 

1-X 2 



+ 2 



X X 



1-X 5 



+ 2e 5 



XX 

1-X 2 



+ 3ee 



/X-X\ 3e/X-X" 



Vi -x 2 



e V 1 - X 2 



X' 



As in the conformal gauge, the expressions for f self and /°' self ? facilitate the well 
renormalisation of the equations of motion (35) and coordinate energy density, 



(38) 
■known 



MA) X- 



1 /X' 



?rad 



MA)e = / 



O.rad 



(39) 
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where the expressions for f rad and f°' rad represent the finite radiation backreaction force. 
The renormalised versions of (36) are 



(40) 



Differentiating (40), gives the power and force due to radiation backreaction 



E = 



MA) 



doe = 



daf 



O.rad 







eX + eX 


-\s> 





^jrad _|_ y0,rad-j£ 



(41) 



Again, some of the terms in (37) and (38) can be shown to be sub-dominant, as for 
the conformal gauge. In particular it is possible to approximate HJj^p, f rad and f°> Tad 
by 



/ 



3X 4 

frad _ 47T/ a 2 A 
3 

0,rad _ 47T/ a 2 A 



eX-ir^lX' 



e V 1 - X 5 



(42) 



e 2 X-X 
1 -X 2 



Substituting the expressions for f°' rad into the power expression (41) yields 

dP 



dl 



3L 



1-X 2 



(43) 



2.7 Eliminating 'runaway' solutions 

This simplified form of the equations of motion using (42) still has serious shortcomings 
because of the presence of the X term. The equations have unphysical, exponentially 
growing or 'runaway' solutions which will, for example, plague any potential numerical 
applications. Furthermore, one would be required to store information at three different 
timesteps, fundamentally changing the nature of a numerical algorithm. It appears, 
however, that both these problems can be circumvented by resubstituting the equations 
of motion, that is, we make the approximations X ps e~ 1 (X'/e) / and X ss e~ 1 (X'/e) / 
in (42) (note we have used the unperturbated equations with e w 0). The equations of 
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Figure 3: The relative positions of the curves y = f(m) (solid line) and y = g(m) (dotted 
line) for typical values of the parameters £1 and OL. Notice that the real positive solutions of 
f(m)=0 — corresponding to the exponentially growing solution of the equations of motion — is not 
a solution of g(m) = 0. 

motion then acquire an analogue of a viscosity term for which there are only damped 
solutions. After performing this resubstitution one finds that the approximate force (42) 
becomes 



/ 



,rad „ 47T/ a 2 A 



0,rad _ 47r/ a 2 Ae 




1 f X.'- (X'/e)' 
1 -X 2 



X.(XVe)' 



(44) 



1 -X 2 

The reason for the suppression of the exponentially growing solution becomes ap- 
parent if we consider simplified one-dimensional equations, 



X - X" = aX , 



X - X" w aX" , 



(45) 



where we have performed the resubstitution assuming that a is small. We now take an 
approximately periodic solution, X" ps — 2 X, and we substitute the ansatz X ~ e mt . 
The solutions for (45) are given respectively by the roots of the following polynomials in 

m, 

f(m) = am 3 - m 2 - O 2 , — ► g(m) = -m 2 - atfm - O 2 . (46) 
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Figure 4: A schematic of the contributions to the radiation backreaction force for a perturbed 
long string configuration, (a) The force is calculated by integrating along O", that is, summing 
the contributions of all string segments on the backward light cone, (b) The expected appearance 
of the actual radiation force contributions; this has a local maximum and a fairly rapid fall-off 
with cancellations. Here, L is the typical wavelength of perturbations on the string (or the string 
curvature radius) and the area under the curve is the total magnitude of the force, (c) The radiation 
force is estimated in the local backreaction approximation using the local magnitude of the force 
and an effective width A for which the area under the two curves is equal. For this approximation 
to be valid, A must be less than the string curvature radius. 

If we rewrite /(to) = (to 2 + Am + Q 2 +B)(am-C), then we see that A = afi 2 + C(a 2 ), 
B = 0{a 2 ) and C = l + 0(a 2 ). If we ignore terms C(ct 2 ), then the solutions of g(m) = 
are approximately solutions of /(to) = 0. However, the real positive solution of /(to) = 0, 
corresponding to the exponentially growing solution of the equations of motion, is not 
a solution of g(m) = 0. Fig. 3 shows the relative positions of the curves y = /(to) and 
y = g{m) for typical values of the parameters Q and a. 

2.8 Understanding the 'local backreaction approximation' 

The local backreaction approximation effectively reduces the problem of calculating the 
backreaction force from the string to the equivalent problem for the electron. The as- 
sumptions underlying (32) and (39) are that the dominant contributions to the integrals 
(27), (28) and (30) come from string segments close to the point under consideration. 
In fig. 4 we have schematically illustrated the construction of the radiation backreac- 
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tion force in the local backreaction approximation. To calculate this force at the time 
t at a particular point on the string (say a = 0), we must sum over the retarded time 
contributions from all other string segments. For definiteness, let us suppose we are con- 
sidering a long straight string with perturbations of typical wavelength L (comparable 
to the string curvature radius R). As we integrate along the backward light cone of 
fig. 4a, we can expect force contributions to take a form appearing something like fig. 4b. 
The precise rate of the fall-off of this force density away from a = is unknown, but 
finiteness certainly implies that it is faster than 1/cr. Moreover, regions of the string 
beyond the curvature radius R will give negative, as well as positive, contributions and 
the resulting net cancellations should ensure rapid convergence of the integrated force at 
large distances. 

The total area under the curve in fig. 4b represents the exact magnitude of the 
radiation backreaction force. The 'local backreaction approximation' to this force is 
illustrated in fig. 4c. We calculate the actual magnitude of the force at the point in 
question and we assume that the contribution from neighbouring segments falls away 
rapidly beyond an effective width A. We then normalize A to ensure that the area 
under the curve (c) is equal to that under curve (b). In §3 we shall discuss the procedure 
for achieving this by comparing with analytic and numerical radiation calculations. 

Of course, we do not expect our local force to evolve every string trajectory com- 
pletely accurately. We are, after all, assuming a uniform A, though it is possible to 
improve this first approximation. Furthermore, the force given by (42), which we shall 
use in the numerical simulations of §3, breaks worldsheet covariance and so we should 
anticipate difficulties describing some special 'null' string trajectories. This is because 
we have taken time as a preferred direction in the derivation of the Lienard-Wiechart po- 
tentials. Effectively, we calculate the backreaction force by summing up all contributions 
from points on the string inside the region D given by 



which corresponds to the shaded diamond in fig. 3a. However, there exists a family 




(47) 
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(48) 

T = T R 



of elongated rectangular regions of equal area related to D by Lorentz boosts. If the 
string trajectory has a typical wavelength (or periodicity) then it seems likely that the 
contributions from the two regions will be similar and the inaccuracies in the force should 
cancel out. We anticipate, therefore, that errors in the local backreaction approximation 
should be small for generic long string trajectories and for closed loops where this 'pseudo- 
periodicity' is likely to be evident. 

In essence the validity of the local backrection approximation hinges on whether the 
following integral is small: 

" ^ 2 J|„_,| >A/2 A-Xdf\\A.X\J 

To summarize, we believe that (48) can be neglected for strings in a realistic network 

because the natural long distance force fall-off will be augmented by strong cancellations 

from a random superposition of distant modes. In any case, at the very least, this 

approximation should work in an some 'averaged' sense. 

2.9 Analytic models for radiative decay 

We can begin to develop confidence in the veracity of the local backreaction approx- 
imation by demonstrating that it predicts the correct scale-dependence of the overall 
radiation power from closed loops and long string trajectories. By analogy with the 
simple model for backreaction of ref. [5], we shall deduce expressions for the evolution of 
the invariant string length L for loops and the relative amplitude £ for perturbed long 
strings, using some of the solutions we presented in §2.2. 

(i) Closed loop solutions 

For closed loops of invariant length L, we can estimate that X ~ 0(1) and X ~ 0(L~ 2 ) 
(the actual average over one period is (|X|) = l/y/2). If we take A ~ L in our local 
approximation, then the power per unit length (43) is proportional to which im- 
mediately recovers the well known result that the power loss from a loop is independent 
of its size L [14] . In general, one can write the radiation power as 

p = r a / a 2 = ^, (49) 
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where T a is some factor dependent on the particular loop trajectory, but not its size and 
k is the radiation backreaction scale, assumed to be independent of timef. The radiation 
damping can be modelled by considering the following equations, 



E = /i I dae = \iL , 
Jo 



(50) 



Integrating these equations one obtains, 

L = L -K{t-t ). (51) 

(ii) Long string solutions 

The generic result for a long string solution, parametrized by its wavelength L and 
relative amplitude E, is X ~ 0{£) and X ~ 0(EL~ 2 ). If A ~ L, then the power per 
unit length in (43) is 

where (3 ~ quantifies the overall strength of the radiation. By analogy to the closed 
loops the radiation damping can be modelled by, 

f L 

E = — dae = fj, + a^iE 2 , 

L Jo (53) 
dP _ _dE _ (3£ 2 

dl dt L 

where a is the (order unity) solution-dependent coefficent of £ 2 in the power series 
expansion of 



lJ 



L dX 3 



da 



(54) 



L J da 

The power loss (52) will lead to an exponential decay of the amplitude and oscillation 
energy per unit length E, 

£ = £oexp(-^y, £ = ^ + a^ 2 exp(--|y . (55) 



f Problems involving global strings become intractable if the logarithmic time dependence of the string 
tension is included. In cosmological problems one finds that the logarithm changes by only about one order 
of magnitude over the enormous timescale between string formation and the present day. We will make this 
assumption in all the following calculations. 
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Exponential decay has already been shown to be generic [5] for a realistic situa- 
tion where the left- and right-moving amplitudes are not precisely the same as in (15). 
However, for some of the periodic solutions of §2.2 such as the 'standing wave' solu- 
tions (14) and (17), the amplitude fall-off was shown analytically to be a power law, a 
fact also confirmed in numerical field theory simulations [5]. How can we reconcile this 
apparent contradiction? The answer lies in noting that the local backreaction approx- 
imation only applies to string configurations in which long-range field correlations are 
suppressed beyond the string curvature radius. Thus when we compare the effects of 
the radiation backreaction force with numerical field theory simulations, we must ensure 
that global correlations are suppressed, eliminating artificial situations with large-scale 
periodic coherence. 

3 Numerical Comparisons 

3.1 Numerical methods 

In order to solve the modified equations of motion (39) and (44) numerically, one must 
recast the simplified resubstituted equations of motion into a first-order form accessible 
to numerical solution. Defining a = X' — eX and = X' + eX, the equations of motion 
can be rewritten as 



Using the above, we were able to evolve string trajectories by modifying a total 
variation non-increasing (TVNI) algorithm [22,9] which has already been well-tested for 
string network evolution in an expanding universe. This method relies on the fact that 
the first order equations of motion (56) are in conservative form, if the backreaction force 
is zero. Artificial compression methods are used to prevent the numerical dissipation of 
kinks. Typically, the algorithm maintains the perpendicularity condition X ■ X' = and 




(56) 
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conserves energy to within a few percent over many timesteps. It should be noted that 
the addition of a small backreaction forcing term does not seem to affect the stability 
of the numerical scheme. However, if the backreaction force becomes larger than the 
tension force, then the equations of motion become qualitatively different behaving like 
a diffusion equation rather than a hyperbolic wave equation. The characteristic Courant 
condition for a diffusion equation is very much more restrictive and so stability problems 
will emerge in this regime. One can address this numerical problem by artifically prevent- 
ing the force from becoming too large, that is, the usual procedure of 'force softening'. 
There is a further technical numerical problem because the coordinate energy e decays 
more rapidly at certain points (for example, at the cusps of the Kibble-Turok loops). 
Eventually, this imposes an unacceptably small timestep on the simulation because we 
always require At < eAa everywhere. However, this problem can be solved by a number 
of approaches, including multiple time-stepping in small e regions, by reparametrising 
the string to redistribute e more evenly, or by eliminating such regions through 'point- 
joining' techniques. 

In all the numerical simulations using the radiation backreaction force in this paper, 
we have employed the constant damping coefficent given by 

^-=^0.001L. (57) 
3MA) 

This choice of damping coefficient corresponds to the cosmologically interesting pareme- 
ter range //(A) ~ 100/f with our numerically determined normalization A ~ 0.1L which 
we shall discuss in the next section. When we compare the results of these Nambu string 
simulations with those using the underlying field theory for which fi(A) ~ 5/f , we have 
had to perform a single global rescaling of the time axis in order to take into account 
the different radiation strengths. 

As well as this one-dimensional effective model, we have developed sophisticated 
numerical algorithms to dynamically simulate string configurations in the Goldstone 
model (2) [5,23,24]. We discretize space on a three-dimensional grid with dimensions 
Ni,N2,Ns in the x, y, z directions respectively, solving the rescaled (f a — > 1, A — > 2) 
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Euler-Lagrange equation, 



a 2 $ a 2 $ a 2 $ d 2 $ 



+ $($$ - 1) = 0. 



(58) 



dt 2 dx 2 dy 2 dz 2 



We employ a second-order leapfrog algorithm for the time derivative and fourth-order 
finite difference approximations for the spatial derivatives. In problems where radiation 
is incident on the boundaries, its sensible to use absorbing boundary conditions [5,24]. A 
second order wave equation which annihilates the reflected wave at the x = boundary 
is 



The efficacy of these methods is discussed in some detail in ref. [5]. 

We use a cylindrically symmetric string ansatz to create initial conditions for both 
long string and loop solutions as in ref. [5]. However, a naive application of this ansatz 
artificially creates long-range correlations which do not conform with the assumptions 
underlying the 'local backreaction approximation'. As we have emphasised, general 
configurations that occur in realistic string networks will not have field correlations 
beyond the average curvature radius of the string because of reconnection processes 
and causality constraints. Consequently, we have modified our ansatz for long string 
configurations by numerically suppressing the initial perturbations with the gaussian, 



where R is the string curvature radius, x is the position in question and x s is the nearest 
long string segment. At large distances r»R, therefore, the string fields will approach 
those for a straight string, as we would expect in a general physical context for random 
small-scale structure. 

To make this distinction plain, fig. 5 illustrates the effect of large distance corre- 
lation suppression on the decay of a periodic sinusoidal solution. The suppressed case 
(fig. 5b) can be seen initially to decay more rapidly than the configuration in which 
perturbations in the fields are correlated out to the simulation boundary (fig. 5a); the 




(59) 
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(60) 
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Figure 5: A comparison of the decay amplitude for a sinusoidal solution in full field theory (a) 
without suppression and (b) with suppression of the field at the curvature radius. Notice that the 
initial decay rate of the suppressed configuration is much faster (exponential) than that for the 
unsuppressed configuration (power law) due to the long range correlations of the latter. 

former is exponential decay, while the latter is power law. However, given the periodic 
boundary conditions, this difference does not persist indefinitely because the suppressed 
configuration will causally relax to an unsuppressed one, as the long range fields become 
correlated on larger and larger scales. Given this limitation imposed by the numerical 
grid size, we can only expect to normalize the local backreaction approximation using 
relatively short simulations (or by using non-periodic configurations). 

3.2 Long string configurations 

We have extensively tested the 'local backreaction approximation', using the modified 
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Figure 6: Log-linear plots for the decay of amplitude for (a) a sinusoidal solution, (b) a pure 
left moving helicoidal perturbation, and (c) a helicoidal perturbation with unequal left- and right- 
moving amplitudes. The straight line typifies the exponential decay. 

Nambu equations of motion (56) and the approximate force derived from (42) by directly 
comparing with field theory simulations of radiating strings in the Goldstone model. 

Using our modified Nambu string simulations, we find that exponential decay is 
generic for all the long string configurations discussed in §2. This decay is illustrated 
in fig. 6 for the sinusoidal solution (17), the pure right-moving helicoidal solution (16) 
and the helicoidal solution with unequal left- and right-moving amplitudes (15). Fig. 7 
illustrates the excellent quantitative agreement with the full field theory simulations by 
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Figure 7: Decay of S using the radiative backreaction force (clotted line) and numerical field 
theory simulations (solid line) for (a) a sinusoidal perturbation, (b) a helicoidal perturbation with 
unequal left- and right-moving amplitudes, and (c) a pure left moving helicoidal perturbation. Note 
the excellent quantitative agreement for all three cases. 

direct comparison with the same three long string solutions. Note that these curves 
have not been matched separately; the same backreaction damping coefficient applies for 
each and there has been only a single global rescaling. The agreement persists for the 
longest time for the (generic) unequal left- and right-moving configuration (15) because 
exponential decay is predicted in this case even after field correlations have relaxed at 
large distances. By comparing with the simple backreaction model for exponential decay 
(55), one can use the numerical field theory results to normalize A, that is, we estimate 



A « (0.1±0.02)L 
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for the long string solutions investigated. This is a result for which there are considerable 
uncertainties at this stage, mainly because of the imprecision inherent in our small- 
scale field theory simulations. We had anticipated that A should be normalized to a 
distance near the string radius of curvature R, which for a sinusoidal perturbation is 
R ~ L/4. The fact that A is smaller than R validates the linearized expansion on which 
the 'local backreaction approximation' is based. The normalisation of the value of A 
above can become ambiguous in certain physical contexts, such as a solution with a 
number of different Fourier modes. In this case, we must make a further approximation 
by normalising to the lengthscale which is radiatively dominant. 

We have also applied these numerical approaches to study the kink solution (19) 
of §2. Fig. 8 compares the evolution of a sharp kink in both the local backreaction 
approximation and in a field theory simulation. Backreaction leads to a substantial 
rounding of the kink, in agreement with the intuitive picture described in refs. [25,5]. 
The results are almost indistinguishable except for the computational advantages of the 
former which, in this case, saved a factor of 10 2 in cpu time and 10 4 in allocated memory. 

We have also performed spectral analysis of the modes on the string using techniques 
similar to those used in three dimensions in ref. [5]. The kink itself can be written as 
an infinite series of odd Fourier modes, while the anticipated endpoint, a sinusoidal 
solution, is just a single Fourier mode. Fig. 9 illustrates the mode decomposition of the 
kink solution initially and then at late times. Radiation backreaction causes decay in 
all modes, but the higher harmonics are clearly damped much more strongly, leading to 
the kink 'rounding'. These spectra can be compared to the kink radiation fields shown 
in ref. [5] which demonstrate the same trend. 

3.3 Closed loop solutions 

We have also applied the local backreaction approximation to the study of loop solutions, 
such as the Kibble- Turok loops described in §2.2. In this case, A ~ L is not independent 
of time because the loop shrinks as it decays*. This problem can be circumvented in the 

* For the long string solutions, periodic boundary conditions forced the solution to have a fixed time- 
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Figure 8: Decay of a kink perturbation (£q = 0.9) using (a) the radiation backreaction force 
and (b) numerical field theory. Notice the visible rounding of the kink in both cases. 

case of loops by choosing the cut-off scale A equal to a constant multiple of the total 
invariant string length, 

L= dae, (62) 
Jo 

which is easily calculable within the evolution algorithm described earlier. 

The damped evolution of the special kinky loop solution (12) is shown after several 
oscillations in fig. 10. As the loop shrinks in size, there is discernible 'rounding' due to 
the radiative damping, though it is less pronounced than in the long string kink decay. 
Unfortunately, evolution for this and other loops could not be continued indefinitely be- 
cause a numerical solution to the Courant violation problem in small e regions has yet to 
be implemented (refer to §3.1). However, the observed 'rounding' is at least qualitatively 
in agreement with a previous attempt to study gravitational backreaction in ref. [26] . In 

independent wavelength A ~ L. 
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Figure 9: The time evolution of the Fourier modes of an initial kink configuration (a). Note the 
damping of higher modes after 16 oscillations in (b). 



Figure 10: The time evolution of a kinky loop solution shown initially and after 5 and 10 
oscillations. Note the decrease in loop size and discernible kink 'rounding'. 
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Figure 11: The effect of radiative damping on loop energy for Ct = Kibble-Turok loops with 
(j) = 71"/ 12 (solid line), 7r/3, (dotted line) and 57T/12 (dashed line). Note the expected linear 
decay of the loop length. 

this non-local approach, all the retarded time radiation contributions were accumulated 
for an unperturbed loop trajectory and then these 'corrections' were applied at the end 
of each oscillation period. Unlike the local backreaction approximation, there is little 
prospect of such 'exact' approaches being implemented in network simulations because 
the 0(N 2 ) algorithms require a supercomputer to evolve a single loop. Nevertheless, we 
anticipate future quantitative comparisons with such methods to determine the accuracy 
of our approach. 

The evolution of the energy of some Kibble-Turok loops is illustrated in fig. 11. 
One can readily observe the linear decay of these solutions, as expected from our simple 
backreaction model (51). Notice, however, the oscillatory nature of the decay due to 
stronger radiation when the loop trajectories becomes more convoluted and when cusps 
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appear. It is interesting to note that a preliminary analysis indicates that, while cusp 
velocities are curtailed by backreaction, their periodic reappearance in these particular 
solutions is not actually prevented. 

The overall decay rate in fig. 11 is parameter dependent; the slope yields the back- 
reaction scale k (or T a ) in (49) which is appropriate for the particular loop trajectory. 
Fig. 12 illustrates the (^-dependence of k for the a = loop solutions (11). This is 
qualitatively similar to analytic estimates of the radiation from these loops, though at 
this stage we can only compare to results for gravitational radiation [14]. Note, how- 
ever, that the divergences at small and large <fi become weakened relative to the previous 
analysis; this may reflect a shortcoming of our approximation or the genuine influence of 
backreaction. However, a considerably more detailed quantitative analysis is necessary 
to test the accuracy of this approach, especially if we are to normalize it properly for 
string network simulations. 

A study of the overall Kibble-Turok loop parameter space (which previous anal- 
yses have regarded as fairly typical) yielded an approximate value k ~ 0.1, given the 
assumed damping coefficient (57) which was set by normalizing A with the long string 
results. However, k 0.1 is the typical backreaction scale expected for GUT-scale global 
string loops with //(A) ps 100/„ (refer to ref. [14,10]), thus independently validating our 
previous normalization 

AwO.lL. (63) 

Given that the case for the local backreaction approximation is not as clear-cut for closed 
loop solutions, these results must be regarded as encouraging. At the very least, this 
approach can be used to phenomenologically incorporate expected loop decay rates, but 
results to date suggest it will do substantially better. 

4 Radiative backreaction in an expanding universe 

Our current understanding of the evolution of a cosmic string network is based on a mar- 
riage between analytic models and sophisticated network simulations [8,9]. However, the 
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Figure 12: Parameter dependence of the backreaction scale K for OL = Kibble-Turok loops 
using the local backreaction approximation. 

network simulations only evolve the free equations of motion for a string in an expanding 
universe. In order to incorporate the radiative effects discussed in the preceding section 
one must modify the equations of motion to include a radiation damping term, 



X + ^(1_X 2 )X- -(— 
a e V e 



= f, 



• , 2h Y2 

e H eX 

a 



(64) 



where a is the scale factor. 

However, to calculate this radiation damping term one must use Green functions in 
an expanding background. In the radiation-dominated era, the retarded Green function 
is 



D ret ( X ,x>) = ^_5((x-xr)e(v-v'). 



(65) 



33 



where x M = (rj, x) and x /fl = (77', x'). In the matter-dominated era the retarded Green 
function also includes the effects of back-scatter off the background spacetime curvature, 

a{rj) 



D ret (x,x') 



S((x - x'f)0(v ~ V) + 1^0(V ~ V' ~ |x - x'| 



(66) 



2ira(r]') 

Applying, the 'local backreaction approximation' in either of these scenarios, one finds 
that the forcing terms are given by, 

f rad = fflat +-g + 0(lA 3 ), 

ct 



/°' rad = / fl °at ad + V + ^(lA 3 ) , 
Ct 



where the flat suffix denotes the flat space backreaction force given by (42) and (g , g) 
is a correction to the force due to the expanding background. Notice that to eliminate 
'runaway' solutions due to X , we must now resubstitute the damped expanding universe 
equations of motion (64). Recall that perturbations with lengthscales r»i/ _1 are essen- 
tially 'frozen' by Hubble damping. It is for the same reason that large-scale perturbations 
will not radiate, despite the high degree of string curvature on these lengthscales. 

The forced equations of motion (64) can be averaged to derive equations for the 
evolution of the density of long strings (poo) and loops (pl), under the influence of the 
expansion, Hubble damping and the radiation backreaction force. If one now inserts a 
term to take into account of loop production, the equations become 



2d 2 cp^ dp 
(1 + (v ))p 00 - 



00 



n a (67) 
Sa cpoo 

PL = PL + , 

a L 

where 

d = d + -d 1 + O(l/t 2 ), 

(X 

(v 2 ) is the average string velocity, d ,di,.. are constants and c is a measure of the 
efficiency of loop production. Substituting p^ = p(/t 2 and L = C,~ l ^ 2 t into (67), one 
obtains 

C 1 



C t 



2-2P(l + (v 2 ))-(c + d)C 1/2 
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where (3 is determined by the scalefactor, a oc t 13 . Eqn (68) has an attractive fixed point, 
which corresponds to the scaling regime. If di = for i > 0, then we have 

c = C 1/2 (l-(v 2 ))-d . 

In the case where d\ is non-zero, one should observe transient effects in the scaling. 
However, for large times these effects will be negligible and the attractive fixed point is 
exactly that for d± = 0. 

5 Conclusions 

We have introduced, and we have endeavoured to justify, a new approach to the study 
of radiation backreaction on strings (and other extended objects and membranes). If 
our analysis and the supporting evidence is valid, then the 'local backreaction approxi- 
mation' offers the hope of quantitative insight into the essentially intractable problem of 
radiative damping effects during string network evolution. It is appropriate, therefore, 
to summarize the main points in our discussion. 

By exploiting the analogy with classical electrodynamics we have used Green func- 
tion methods to separate the self- and radiation fields of a global string, using the former 
to renormalize the string tension. We then approximated the radiation force at a point 
on the string by an expansion in powers of a cut-off parameter A. The 'local back- 
reaction approximation' to the radiation force, then, is the local force at the point in 
question multiplied by an effective width A beyond which the neighbouring string seg- 
ment contributions become negligible. We normalize A in order to reproduce the actual 
radiation force in known situations and we have confirmed the self-consistency of the 
approximation by demonstrating that A is less than the string curvature radius. We 
then generalized this approach to the temporal transverse gauge, a convenient gauge for 
studying dissipative string processes. The final step was to remove unphysical 'runaway' 
solutions, which plague the analogous point-particle analysis, by resubstituting the string 
equations of motion. 
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We then tested a numerical implementation of the 'local backreaction approxima- 
tion' by investigating a variety of long string and closed loop trajectories. We have 
directly compared the results with analytic radiation calculations and numerical field 
theory simulations, demonstrating a consistent normalization for A using these inde- 
pendent methods. This approach reproduces the correct scale-dependence of radiative 
effects and demonstrates satisfactory quantitative agreement for a wide variety of dif- 
ferent solutions. There is clearly scope for a more detailed analysis of the accuracy of 
this approach and for addressing a number of outstanding issues. However, we have 
presented sufficient grounds for believing that the 'local backreaction approximation' is 
a significant step forward in the study of string radiative backreaction. 
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